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We study properties of topological phases by calculating the ground state degeneracy (GSD) of 
the 2d Levin- Wen (LW) model. Here it is explicitly shown that the GSD depends only on the spatial 
topology of the system. Then we show that the ground state on a sphere is always non-degenerate. 
Moreover, we study an example associated with a quantum group, and show that the GSD on a torus 
agrees with that of the doubled Chern-Simons theory, consistent with the conjectured equivalence 
between the LW model associated with a quantum group and the doubled Chern-Simons theory. 
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I. INTRODUCTION 

In recent years two-dimensional topological phases 
have received growing attention from the science com- 
munity. They represent a novel class of quantum mat- 
ter at zero temperature whose bulk properties are ro- 
bust against weak interactions and disorders. Topological 
phases may be divided into two families: doubled (with 
time-reversal symmetry, or TRS, preserved), and chiral 
( with TRS broken). Either type may be exploited to do 
fault-tolerant (or topological) quantum computing^—. 

Chiral topological phases were first discovered in inte- 
ger and fractional quantum Hall (IQH and FQH) liquids. 
Mathematically, their effective low-energy description is 
given by Chern-Simons theorj^ or (more generally) topo- 
logical quantum field theory (TQFT)£. One character- 
istic property of FQH states is ground state degeneracy 
(GSD), which depends only on the spatial topology of the 
systemZr— and is closely related to fractionizationiS - — 
of quasiparticlc quantum numbers, including fractional 
(braiding) statistic a 13 ' 14 . In some cases the GSD has been 
computed in refs i 12 ' 15 . 

Chern-Simons theories are formulated in the contin- 
uum and have no lattice counterpart. Doubled topolog- 
ical phases, on the other hand, do admit a discrete de- 
scription. The first known example was Kitaev's toric 
code models. 

More recently, Levin and Wen (LW) 1 ^ constructed 
a discrete model to describe a large class of doubled 
phases. Their original motivation was to generate 
ground states that exhibit the phenomenon of string-net 
condensation^ as a physical mechanism for topological 
phases. The LW model is defined on a trivalent lattice (or 
graph) with an exactly soluble Hamiltonian. The ground 
states in this model can be viewed as the fixed-point 
states of some renormalization group flo w 18 ' 19 . These 
fixed-point states look the same at all length scales and 
have no local degrees of freedom. 

The LW model is believed to be a Hamiltonian ver- 
sion of the Turaev-Viro topological quantum field the- 
ory (TQFT) in three dimensional spacetime 4 - ' 20 ' 21 and, 
in particular cases, discrctized version of doubled Chern- 



Simons theor y 22 ' 23 . Like Kitaev's toric code model^, we 
expect that the subspace of degenerate ground states in 
the LW model can be used as a fault-tolerant code for 
quantum computation. 

In this paper we report the results of a recent study on 
the GSD of the LW model formulated on a (discrctized) 
closed oriented surface M. Usually the GSD is exam- 
ined as a topological invarian t 20 ' 21 ' 23 of the 3-manifold 
S 1 x M. In a Hamiltonian approach accessible to physi- 
cists, we will explicitly demonstrate that the GSD in the 
LW model depends only on the topology of M on which 
the system lives and, therefore, is a topological invari- 
ant of the surface M. We also show that the ground 
state of any LW Hamiltonian on a sphere is always non- 
degenerate. Moreover, we examine the LW model as- 
sociated with quantum group Sf7fc(2), which is conjec- 
tured to be equivalent to the doubled Chern-Simons the- 
ory with gauge group SU (2) at level k, and compute the 
GSD on a torus. Indeed we find an agreement with that 
in the corresponding doubled Chern-Simons theory^21. 
This supports the above-mentioned conjectured equiva- 
lence between the doubled Chern-Simons theory and the 
LW model, at least in this particular case. 

The paper is organized as follows. In Section II we 
present the basics of the LW model, easy to read for 
newcomers. In Section III topological properties of the 
ground states are studied, and the topological invariance 
of their degeneracy is shown explicitly. In section IV we 
demonstrate how to calculate the GSD in a general way. 
In section V we provide examples for the calculation par- 
ticularly on a torus. Section VI is devoted to summary 
and discussions. The detailed computation of the GSD 
is presented in the appendices. 



II. THE LEVIN- WEN MODEL 

Start with a fixed (connected and directed) trivalent 
graph r which discretizes a closed oriented surface M 
(such as a torus). To each edge in the graph we assign a 
string type j, which runs over a finite set j — 0, 1, N. 
Each string type j has a "conjugate" j* that describes 
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the effect of reversing the edge direction. For example j 
may be an irreducible representation of a finite group or 
(more generally) a quantum group^. 

Let us associate to each string type j a quantum dimen- 
sion dj, which is a positive number for the Hamiltonian 
we define later to be hermitian. To each triple of strings 
{i,j,k} we associate a branching rule Sijk that equals 1 
if the triple is "allowed" to meet at a vertex, if not 
(in representation language the tensor product i® j ®k 
cither contains the trivial representation or not). This 
data must satisfy (here D = V ■ $) 



^ dkSijk* = didj 

k 

d.d,!),, i.- = d k D 



(1) 



j = is the unique "trivial" string type, satisfying 0* = 



and 5, 
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The Hilbert space is spanned by all configurations 
of all possible string types j on edges. The Hamilto- 
nian is a sum of some mutually-commuting projectors 
H := — J2 V Qv — Yip Bp (one for each vertex v and each 
plaquette p). Here each projector Q v = S^k with i, j, k on 
the edges incoming to the vertex v. Q v = 1 enforces the 
branching rule on v. Throughout the paper we work on 
the subspace of states in which Q v = 1 for all vertices. 
Each projector B p is a sum D~ 1 Y, s dsBp of operators 
that have matrix elements (on a hexagonal plaquette for 
example) 



(2) 
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Here Vj = \fdj is real. The symmetrized 6j symbol: 
are complex numbers that satisfy 
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pentagon id: 
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For example, these conditions are known to be 
satisfied^ if we take the string types j to be all irre- 
ducible representations of a finite group, dj to be the 
dimension of corresponding representation space, and G 
to be the symmetrized Racah 6j symbols for the group. 
In this case the LW model can be mapped 2 ^ to Ki- 
taev's quantum double model^. More general sets of data 
{G,d, 6} can be derived from quantum groups (or Hopf 
algebras)^. We will discuss such a case later using the 
quantum group 5f7fc(2) (fc being the level). 





p(2) 



FIG. 1: Given any two trivalent graphs T' 1 - 1 and F^ 2 ' dis- 
cretizing the same surface, we can always mutate r' 1 ' to r' 2 ' 
by a composition of elementary / moves. In general F' 1 ' and 
F^ 2 ' are not required to be regular lattices. Fhese diagrams 
happen to be the same as2&, but in a slightly different context. 



III. GROUND STATES 

Any ground state |<f>) (there may be many) must be a 
simultaneous +1 eigenvector for all projectors Q v and B p . 
In this section we demonstrate the topological properties 
of the ground states on a closed surface with non-trivial 
topology. 

Let us begin with any two arbitrary trivalent graphs 
rW and r^ 2 ) discretizing the same surface (e.g., a torus). 
If we compare the LW models based on these two graphs, 
respectively, then immediately we see that the Hilbert 
spaces arc quite different from each other (they have dif- 
ferent sizes in general). 

However, we may mutate between any two given triva- 
lent graphs rW and T^ 2 ) by a composition of the follow- 
ing elementary movea 2 ^ (see also FigQ]): 

=> J , for any edge; 



, for any vertex. 



, for any triangle structure. 



Suppose we are given a sequence of elementary / moves 
that connects two graphs T^ 1 ' — ^ T^ 2 ). We now construct 
a linear transformation H 1 - 1 ^ — s- Ti!^ between the two 
Hilbert spaces. This is defined by associating linear maps 
to each elementary / move: 
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The mutation transformations between T-L^ and T-L^ 
are constructed by a composition of these elementary 
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maps. As a special example, the operator B p = 
D^ 1 J2 s dsBp is such a transformation. In fact, on 
the particular triangle plaquette p as in ([?]), we have 
-B p =v = T2T3, by using the pentagon identity in ([3]). 

Mutation transformations are unitary on the ground 
states. To see this, we only need to check that the el- 
ementary maps Ti, T2, and T 3 are unitary. First note 
that the following relations hold: f\ = f u T 2 f = T 3 , and 
T3 = T2. We emphasize that these are maps between the 
Hilbcrt spaces on two different graphs. For example, we 
check t\ = T\ by comparing matrix elements 
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where in the third equality we used the symmetry condi- 
tion in © • 

Similarly, for Tj = T 3 (or T3 = T2), we have 
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Now we verify unitary. First, T^Ti = id and T 2 T 2 = 
T3T2 = id by the orthogonality condition in ([3]) (note 
that, since we have not used any information about the 
ground states in this argument, T\ and T2 are unitary 
on the entire Hilbert space). For unitary of T 3 we check 



T 3 T T 3 = T 2 T 3 



1. The last equality only holds on the 
ground states since we have already seen that T2T3 = 
B p= y and B p= y = 1 only on the ground states. 

As another consequence of the above relations, the 
Hamiltonian is hcrmitian since all _B p 's consist of ele- 
mentary 7\, T 2 , and T 3 maps. Particularly, on a triangle 



plaquette, we have B p= 



p=V 



(T 2 T 3 )t 



Tin 



T 2 T 3 



The mutation transformations serve as the symmetry 
transformations in the ground states. If |<f>) is a ground 
state then T|$) is also a ground state, where T is a 
composition of Tj's associated with elementary / moves 
from r^ 1 ) to r' 2 - 1 . This is equivalent to the condition 
t (U p b p) = (Up'B p ')f, which can be verified by the 



conditions in Q . (Here p and p' run over the plaquettes 
on and T^ 2 > , respectively. Also note that the -Bp's are 
mutually-commuting projectors, i.e., B p B p = B v , and 
thus Y[ p Bp is the projector that projects onto the ground 
states.) 

These symmetry transformations look a little different 
from the usual ones since they may transform between 
the Hilbcrt spaces TiS 1 ^ 1 and H.^ on two different graphs 
rW and r( 2 ). In general, and do not have the 
same number of vertices and edges. And thus and 
have different sizes. However, if we restrict to the 
ground-state subspaces T-L^ and "Hg 2 \ mutation trans- 
formations are invertible. In fact, they are unitary as we 
have just shown. 

The tensor equations on the 6j symbols in ([3]) give 
rise to a simple result: each mutation that preserves the 
spatial topology of the two graphs induces a unitary sym- 
metry transformation. During the mutations, local struc- 
tures of the graphs are destroyed, while the spatial topol- 
ogy of the graphs is not changed. Correspondingly, the 
local information of the ground states may be lost, while 
the topological feature of the ground states is preserved. 
In fact, any topological feature can be specified by a topo- 
logical observable O that is invariant under all mutation 
transformations f from to H^: O'f = TO (where 
O is defined on the graph and O' on r' 2 )). 

The symmetry transformations provides a way to char- 
acterize the topological phase by a topological observ- 
able. In the next section we will investigate the GSD as 
such an observable. 

Let us end this section by remarking on uniqueness of 
the mutation transformations. There may be many ways 
to mutate to T^ 2 ) using fx, and / 3 moves. Each 
way determines a corresponding transformation between 



(g) the Hilbert spaces of ground states, h[P and Wq'' . It 
turns out that all these transformations are actually the 
same if the initial and final graphs to r*- 2 -* are fixed, 
i.e., independent of which way we choose to mutate the 
graph to r^ 2 ). This means that the ground state 
Hilbert spaces on different graphs can be identified (up 
to a mutation transformation) and all graphs are equally 
good. 

One consequence of the uniqueness of the mutation 
tranformation is that the degrees of freedom in the 
ground states do not depend on the specific structure 
of the graph. In this sense, the LW model is the Hamil- 
tonian version of some discrete TQFT (actually, Turaev- 
Viro type TQFT, see^I). The fact that the degrees of 
freedom of the ground states depend only on the topol- 
ogy of the closed surface M is a typical characteristic of 
topological phases^"— 1 12 ' 15 . 



,(2) 



IV. GROUND STATE DEGENERACY 

In this section we investigate the simplest nontrivial 
topological observable, the GSD. Since JT B p is the pro- 
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jector that projects onto the ground states, taking a trace 
computes GSD = tr(FJ B p ). 

We can show that GSD is a topological invariant. 
Namely, in the previous section we mentioned that, 
by using ([3]), Y[ p Bp is invariant under any muta- 
tion T between the Hilbcrt spaces T-LS 1 ^ and HS 2 ^ : 
fHHp> B P ')f = Up b p- Taking a trace of both sides 
leads to tr'(FJ , B p >) = tr(FJ S p ), where the traces are 
evaluated on HS 2 ^ and H.^ respectively. 

The independence of the GSD on the local structure of 
the graphs provides a practical algorithm for computing 
the GSD, since we may always use the simplest graph 
(see Fig [2] and examples in the next section) . 

Expanding the GSD explicitly in terms of 6j symbols 
using ([2]) we obtain 



GSD 
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(7) 

The formula needs some explanation. P is the total num- 
ber of plaquettes of the graph. Each plaquette p con- 
tributes a summation over s p together with a factor of 

^y- . In the picture in © the top plaquette is being oper- 
ated on first by Bp^ , next the bottom plaquette by B^ , 
third the left plaquette by Bp^ , and finally the right pla- 
quette by Bp* . Although ordering of the Bt, operators is 

not important (since all -Bp's commute with each other), 
it is important to make an ordering choice (for all pla- 
quettes on the graph) once and for all. 

Each edge e contributes a summation over j e and j' e 
together with a factor of dj e dj> . Each vertex contributes 
three 6j symbols. 

The indices on the 6j symbols work as follows: since 
each vertex borders three plaquettes where -Bp's are being 
applied, we pick up a 6j symbol for each corner. However, 
ordering is important: because we have an overall order- 
ing of Bp's, at each vertex we get an induced ordering 
for the 6j symbols. Starting with the 6j symbol furthest 
left we have no primes on the top row. The bottom two 
indices pick up primes. All of these variables (primed or 
not) are fed into the next 6j symbol and the same rule 
applies: the bottom two indices pick up a prime with the 
convention ()" = (). 

By the calculation of the GSD, we have characterized 
a topological property of the phase using local quantities 
living on a graph discrctizing AI of nontrivial topology. 





(a) 



FIG. 2: All trivalent graphs can be reduced to their simplest 
structures by compositions of elementary / moves, (a) on a 
sphere: 2 vertices, 3 edges, and 3 plaquettes. (b) on a torus: 
2 vertices, 3 edges, and 1 plaquette. 



V. EXAMPLES 

(1 ) On a sphere. To calculate the GSD, we need to input 
the data {G z ^,dj, Sij m } and evaluate the trace in ([7]). 
We start by computing the GSD in the simplest case of 
a sphere. 

Let's consider the simplest graph as in Fig. [2ja). 
We show in Appendix A that the ground state is non- 
degenerate on the sphere without referring to any specific 
structure in the model: GSD spherG = 1. In fact, for more 
general graphs one can write down^ the ground state as 
FJ p -Bp 1 0) up to a normalization factor, where in |0) all 
edges are labeled by string type 0. 

We notice that the GSD on the open disk (which is 
topologically the same as the 2d plane) can be studied 
using the same technique. This is because the open disk 
can be obtained by puncturing the sphere in Fig [2](a) at 
the bottom. Although this destroys the bottom plaque- 
tte, we notice that the constraint B p = 1 from the bottom 
plaquette is automatically satisfied as a consequence of 
the same constraint on all other plaquettes. The fact 
that GSD sphoro (= GSD disk ) = 1 indicates the non-chiral 
topological order in the LW model. 

(2) Quantum double model. When the data are deter- 
mined by representations of a finite group G, the LW 
model is mapped to Kitaev's quantum double modeli^. 
The ground states corresponds one-to-one to the flat G- 
conncctionsi. The GSD is 



GSDqd — 



Hom(7n (M),G) 



G 



(8) 



where Hom(7Ti(A^), G) is the space of homomorphisms 
from the fundamental group iti(M) to G, and G in the 
quotient acts on this space by conjugation. 
In particular, the GSD ([8]) on a torus is 



GSD'q ™ 8 



|{(o,6)|o,6G G-^ba^b- 1 = e}/ ~| (9) 

where ~ in the quotient is the equivalence by conjugation, 

(a, b) - {hah- 1 , hbh- 1 ) for all h e G 

The number ([9]) is also the total number of irreducible 
representations^ of the quantum double -D(G) of the 
group G. On the other hand, the quasiparticles in the 
model arc classified^ by the quantum double -D(G). Thus 
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the GSD on a torus is equal to the number of particle 
species in this example. 

(3) SUk(2) structure on a torus. More generally, on a 
torus any trivalent graph can be reduced to the simplest 
one with two vertices and three edges, as in Fig[Ub). On 
this graph the GSD consists of six local 6j symbols. 

GSD " X! ,l ~' l r< 1 ,, d n < l ,\ •'/, •'/. 

/ rihhil f~ii'i 'jij' 3 n hH ii \ ( /-tiZhi* r>Hii\ riili'-z h \ 

(10) 

Now let us take the example using the quantum group 
SUk(2). It is known that SUk(2) has k + 1 irreducible 
representations, and thus the GSD we calculate is finite. 
We take the string types to be these representations, la- 
beled as 0,1,..., k, and the data dj, 6ij m } to be 
determined by these representations (for more details, 
seeM^). 

In Appendix B we show that in this case (for the LW 
model on a torus with string types given by irreps of 
SU k {2)) we have GSD = (k + l) 2 . We argue this both 
analytically and numerically. 

On the other hand, it is widely believed that when the 
string types in the LW model are irreps from a quantum 
group at level k, then the associated TQFT is given by 
doubled Chern-Simons theory associated with the corre- 
sponding Lie group at level zfc fe 24 ' 32 . This equivalence 
tells us that in this case the LW model can be viewed as 
a Hamiltonian realization of the doubled Chern-Simons 
theory on a lattice, and it provides an explicit picture of 
how the LW model describes doubled topological phases. 

Along these lines, our result is consistent 3 - 3 - with the 
result GSDcs = k + 1 for Chern-Simons SU(2) theory 
at level A; on a torus. This can be seen since the Hilbert 
space associated to doubled Chern-Simons should be the 
tensor product of two copies of Chern-Simons theory at 
level ±fc. 



VI. SUMMARY AND DISCUSSIONS 

In this paper, we studied the LW model that describes 
2d topological phases which do not break time-reversal 
symmetry. By examining the 2d (trivalent) graphs with 
same topology which are related to each other by a 
given finite set of operations (Pachner moves), we de- 
veloped techniques to deal with topological properties of 
the ground states. Using them, we have been able to 
show explicitly that the GSD is determined only by the 
topology of the surface the system lives on, which is a typ- 
ical feature of topological phases. We also demonstrated 
how to obtain the GSD from local data in a general way. 
We explicitly showed that the ground state of any LW 
Hamiltonian on a sphere is non-degenerate. Moreover, 
the LW model associated with quantum group SUk(2) 
was studied, and our result for the GSD on a torus is 



consistent with the conjecture that the LW model asso- 
ciated with quantum group is the realization of a doubled 
Chern-Simons theory on a lattice or discrete graph. 

Finally, let us indicate possible extension of the results 
to more general cases. First, more generally in the LW 
model, an extra discrete degree of freedom, labelled by an 
index a, may be put on the vertices. Then the branch- 
ing rule 5?j k , when its value is 1, may carry an extra 
index a. (In representation language this implies that 
given irreducible representations i,j and k, there may be 
multiple inequivalent ways to obtain the trivial represen- 
tation from the tensor product of i j ® k. The index 
a just labels these different ways.) The 6j symbols ac- 
cordingly carry more indices. (For more details sec the 
first Appendix in the original paper— of the LW model.) 
The expression {7J for GSD is expected to be generaliz- 
able to these cases. Secondly, the spatial manifold (e.g. 
a torus) on which the graph is defined may carry non- 
trivial charge, e.g. labelled by ii in the SUk(2) case. This 
corresponds to having a so-called fluxon excitation (of 
type ii) above the original LW ground states. The lowest 
states of this subsector in the LW model coincide with 
the ground states for the Hamiltonian obtained by re- 
placing the plaquette projector B p = D^ 1 J2j djB 3 p with 

Bp = D^ 1 J2j SijB 3 pl where is the modular S- matrix. 
(See Appendix B.) The GSD in this case is computable 
too, but we leave this for a future paper—. 
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Appendix A: GSD = 1 on a sphere 

In appendix, we derive GSD = 1 on a sphere for a gen- 
eral Levin- Wen model, without referring to any specific 
structure of the data {d, 8, G}. All we will use in the 
derivation are the general properties in eq. (fTJ) and eq. 
©• 

The simplest trivalent graph on a sphere has three 
plaquettes and three edges, as illustrated in Fig. H^a). 
Following the standard procedure as in (J7J), the GSD is 
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expanded as 



they are 
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where B pi is acting on the top bubble plaquctte, B p2 
on the bottom bubble plaquette, and B P3 on the rest 
plaquette outside the two bubbles. 

All 6j symbols can be eliminated by using the orthog- 
onality condition in eq. (|3|) three times, 

Ur^r'ji'fs r*j 3 *ji ~ ^ °li3^3 °3i323 3 

EJ /~l333l*32 ^33*3231 _ _}_c r 
" slJ »'ji*ir S *J'U2 ~~ d . l 3'i323 3 * 3l3 2 3 3 * 



't*i2J3 2 _ J. 5 3ih3 3 5 313233* ( A2 ) 
"3i 



and the GSD is a summation in terms of {d, (5}: 
GSD p = — ^ d.j' i d. ] ' :2 dj 3 8j 1 j 2 j*8y i: j 2 y 3 *8 : j 1 ji 2 j^ 



3l32333\3 2 33 



(A3) 



Summing over 7^, j'2, and j'3 using (Q} finally leads to 
GSD sphcrc = 1. 



Appendix B: GSD on a torus for SUk(2) 

Let us consider the example associated with the quan- 
tum group SUk{i) (with the level k an positive integer) 
and calculate the GSD on a torus. 

There are fe + 1 string types, labeled as j = 0, 1, 2, k. 
They are the irreducible representations of SUk(2). The 
quantum dimensions dj arc required to be positive for all 
j, in order that the Hamiltonian is hermitian. Explicitly, 



dj = 



(j+l)7T 

fc+2 



sm 



fc+2 



D = Td 2 



k + 2 



^ "' 3 2 sin 2 - 2 - 

j=0 Zbm fe+2 

The branching rule is <5 rst = 1 if 



(Bl) 



r + s + t is even 

r + s > t, s + t > r,t + r > s 

r + s + t<2k 



(B2) 



and 5 rs t = otherwise. The explicit formula for the 
6j symbol can be found i n 29 ' 30 . However, we do not 
need the detailed data of the 6j symbol in the following 
computation of the GSD. 

Let us start with formula in (fT0|) . and reorder the 6j 
symbols, 



GSD =D- X V d 
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v r v r v n v i3 G{\\ 2 1 3 G jl , h ,i 3 

J2 J3 J 2 33 S* 7 3 7, S* 1, 



=d- 1 y d s ( v h v j3 v, v r G j i j *ii &*?'.?;) 

s y 31 33 3-l 3 3 s*j[*3 3 sf 3 3[ J 



S3l32333i3 2 33 



{v^v^G^G^l 
(v A v 3k v 32 v j3 G^G^ 



(B3) 



where the symmetry condition in ([3|) was used in the 
second equality. 

Let us compare the formula in (|B3[) with that in (|A1[) . 
We set j = j* for all j and drop all stars, since all irre- 
ducible representations of SUk(2) are self-dual. Then we 
find that the summation (|B3p has the same form as the 
trace of D^ 1 J2 S d s B p2 B p3 B pi on the graph on a sphere 
as in ED, 



tr l 



: E< 

313233 



-Y 

D ^ 

S 

L 33 



dsK) 



j}^d s B s p2 B s P3 B s Vi 




(B4) 



where B p is defined on the only plaquette p on the torus 

(see Fig. HJb)), while B* Bi 2 B p3 is defined on the same 
graph on a sphere as in (|A1[) (sec Fig. HJa)). 

The GSD on a torus becomes a trace on a sphere. The 
latter is easer to deal with since the ground state on a 
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sphere is non-degenerate. The counting of ground states 
on a torus turns into a problem dealing with excitations 
on the sphere. 

In the following we evaluate the summation in the rep- 
resentation of elementary excitations, let us introduce a 
new set of operators {hp} by a transformation, 



— nL 

S r0 



(B5) 



Here s rs is a symmetric matrix (referred to as the mod- 
ular ^-matrix for SUk(2)), 



Applying this rule reduces the summation (|B4[) to 

S 

=*4e* e 



S sri S sr2 S sr3 ~r 2 ~r 3 ~r^ 

„ „ „ U P2 H P3 ' l Pl > 
s ri0 s r 2 s r 3 



y - 



d 



S ri 0Sr 2 0Sr 3 



J rir 2 r 3 



(BIO) 



Then we substitute (fBTj), (|B2|) and (JB6J) in and obtain 



■ (r+l)(s+l)ir 
1 Sm fc+2 



sin 



(B6) 



A;+2 



and has the properties 



E 



5 sr , S r = dr/VD 



8rs 



(B7) 



Eq. (|B5[) can be viewed as a finite discrete Fourier 
transformation between {h r p } and {B^}. By properties 
(|B7|> . we see that {n^} are mutually orthonormal projec- 
tors, and they form a resolution of the identity: 



S rs n r p , 



id 



(B8) 



In particular, h p = j)^Z s d s B p is the operator B p in 
the Hamiltonian. The operator hp projects onto the 
states with a quasiparticlc (labeled by r type) occupy- 
ing the plaquctte p. Expressed as common eigenvectors 
of {hp}, the elementary excitations are classified by the 
configuration of these quasiparticles. 

Particularly, on the graph on a sphere as in (|B4|) . the 
Hilbert space has a basis of {|ri, v?,, 1*3)}, where only those 



T\ , r 2 , and r 3 that satisfy 5 r 



1 are allowed. Each 



basis vector |ri, 7*2^3) is an elementary excitation with 
the quasiparticles labeled by n, r2, and r% occupying the 
plaquettes p\ , pi , and p^ . The configuration of quasipar- 
ticles are globally constrained by 5 ri r 2 r a = 

1.34. Thcrc _ 

fore, tracing opertors {h r p } leads to 



QT^torus 
^ U SU k (2) 



k 

E 



(ri + l)7T • (r2 + l)7T . (r3 + l)7T 

ri,r 2) r 3 =0 Sm k+2 Sm fc+2 Sm fc+2 



Sill 



fe+2 



A: r 

EE: (r+l)7r • (s+l)7r • (r — .s+l)7r 
r=0 s=0 Sm fe+2 S111 fc+2 Sm fc+2 



= (fc+l) 



(Bll) 



(Here we omit a rigorous proof of the last equality.) 

We can also verify GSD = (fc+1) 2 by a direct numerical 
computation. We take the approach in 3 ^ to construct the 
numerical data of 6j symbols. The construction depends 
on a parameter, the Kauffman variable A (in the same 
convention as in 3 ^) , which is specialized to roots of unity. 
We make the following choice: 



A = exp(7ri/3) at k = 1 
A = exp(3m/8) at k = 2 
A = exp(37ri/5) at k = 3 



(B12) 



By this choice, the quantum dimensions dj take the 
values as in (|B1[) . and the 6j symbols satisfy the self- 
consistent conditions in ([3J. Using such data of quantum 
dimensions dj and 6j symbols, We compute the summa- 
tion (TlOl) at 



GSD = 4 at k = 1 
GSD = 9 at k = 2 
GSD = 16 at k = 3 



(B13) 



(B9) which verifies GSD = (k + l) 2 in the particular cases. 
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